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Abstract 

Lepage's improvement scheme is a recent major progress in lattice QCD, 
allowing to obtain continuum physics on very coarse lattices. Here we dis- 
cuss improvement in the Hamiltonian formulation, and we derive an improved 
Hamiltonian from a lattice Lagrangian free of 0(a 2 ) errors. We do this by the 
transfer matrix method, but we also show that the alternative via Legendre 
transformation gives identical results. We consider classical improvement, 
tadpole improvement and also the structure of Liischer-Weisz improvement. 
The resulting color-electric energy is an infinite series, which is expected 
to be rapidly convergent. For the purpose of practical calculations, we con- 
struct a simpler improved Hamiltonian, which includes only nearest-neighbor 
interactions. 
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1 Introduction 



Lepage's improvement scheme[||] is a recent major progress in lattice QCD, 
opening the possibility to approach continuum physics on very coarse lattices. 
In this paper we want to address the question how this improvement scheme 
can be formulated for a lattice gauge theory in the Hamiltonian approach. 
Although standard lattice gauge theory has been very successful over the last 
two decades, there are areas where progress has been quite slow. Examples 
are the dynamical computation of the ^-matrix and cross sections, QCD at 
finite baryon density, or the computation of QCD structure functions in the 
region of small xb and Q 2 . This situation calls for the development of new 
methods, and in our opinion the lattice Hamiltonian approach is a viable 
alternative || which should be explored. The Hamiltonian approach corre- 
sponds to consider a continuous time, i.e. at = 0. Similar ideas have been 
pursued recently by workers in standard lattice gauge theory by considering 
anisotropic lattices with lattice spacings a t << a s . This has the purpose to 
improve the computation of the mass spectrum |], |3|. 

As a first step, we here restrict ourselves to the problem of improvement 
in pure gauge theory. Let us recall that the improvement program by Lepage 
consists of several steps, by starting from the Wilson action: First defining 
a classically improved action, second performing tadpole improvement and 
third introducing additional quantum corrections (Luscher-Weisz improve- 
ment). We now discuss how to carry over these ideas to the Hamiltonian 
formulation. Different strategies are possible, let us explain those for the 
case of classical improvement. 

Strategy 1: Construct the classical Hamiltonian corresponding to the clas- 
sical Wilson action. Improve this classical Hamiltonian and quantize this 
Hamiltonian according to the rules of canonical quantization. This yields a 
classically improved quantum Hamiltonian. 

Strategy 2: Starting from the classically improved Wilson action, construct 
via the transfer matrix a classically improved quantum Hamiltonian. 
Strategy 3: Starting from the classical Wilson action, construct first the cor- 
responding quantum Hamiltonian via the transfer matrix. This yields the 
Kogut-Susskind Hamiltonian H(E, U) where U and E are the link variables 
and their canonical conjugate momenta. The usual expansion in powers 
of the field variables A and their conjugate d/dA yields then the standard 
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expression / d 3 x((d/dA) 2 + B 2 (A)) up to correction of (a 2 ). This Kogut- 
Susskind Hamiltonian may be improved by adding corrections such that an 
better agreement with this formal continuum limit operator is obtained. 

Different strategies exist also for the further improvements (tadpole and 
Luscher-Weisz) with respect to the quantum Hamiltonian. In principle, the 
construction of these quantum corrections should start from new perturba- 
tive calculations in the Hamiltonian framework. The coefficients of a fully 
improved action as given by Lepage can only be used as starting point for 
the transfer matrix to obtain a fully improved Hamiltonian if the action is 
expressed on a lattice with the time spacing being much smaller than the 
spatial spacing (a t << a s ). 

In this paper we restrict ourselves mostly to a discussion of classical im- 
provement. As a first result we show explicitly that the first strategy - the 
canonical quantization of a classical lattice gauge theory - is a viable alterna- 
tive to the second strategy - using the transfer matrix - leading to the same 
quantum Hamiltonian, but in a more direct way. 

The classically improved quantum Hamiltonian obtained in this way has 
the mathematical structure of a kinetic part with an infinite number of terms. 
The reason for this structure is given by the fact that the inversion of a nearly 
local matrix leads to a non-local matrix. This being an undesirable feature 
from the point of view of practical calculations, we show that it is possible 
to use the non-uniqueness of the improved action to obtain an improved 
quantum Hamiltonian containing only nearest neighbor interaction terms. 

Finally we discuss the structure of the quantum Hamiltonian related to 
the tadpole and Luscher-Weisz corrections. The determination of the cor- 
responding coefficients, however, will be deferred to a future investigation. 
Although we do not discuss the third strategy in detail, a fully consistent 
computation of the quantum corrections should eventually lead to the same 
result as the first two strategies. Note that we discuss here only the im- 
proved Hamiltonian for the purpose to compute the spectrum. The general 
existence of such improved Hamiltonians is discussed in standard many-body 
theory in the context of model space calculations ([[l^, Like in the action 
formulation general observables require particular improved operators. 
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Calculations of the glueball spectrum using the coupled cluster method 
based on the standard Kogut-Susskind Hamiltonian have been done by Luo 
et al. H and Schiitte et al. ||. An incorporation of an improved Hamiltonian 
should be possible and one would expect reliable results already in lower order 
of the coupled cluster truncation compared to the standard Kogut-Susskind 



Hamiltonian. Calculations in this direction are under way [10 



2 From Wilson action to Kogut-Susskind 
Hamiltonian 

2.1 Canonical method via Legendre transformation 

Before deriving the improved Hamiltonian, we describe in a pedagogical man- 
ner how to obtain the standard Kogut-Susskind lattice Hamiltonian M from 
the classical lattice Lagrangian using the Legendre transformation [6|, 0], 
and canonical quantization. Wilson's Euclidean lattice action is given by 

s E =°^£(i-Po )+ ^i:(i-Po). (i) 

at 9 2 fcu a 9 t^u 

Here the notation of Ref.|13| has been used, i.e., t — □ stands for time-like 



and s — □ for space-like plaquettes, respectively, and 

P D = ±-ReTr(U D ). (2) 

For later use we need to distinguish between the Euclidean and Minkowski ac- 
tion as well as Lagrangian. Its relation is defined, when going from Minkowski 
to Euclidean time by the transformation it — > t, exp^S^,/] — > exp[— Sg], 
L M — > —L E and Sm = J dtL M as well as Se = J dtL E . Thus the Euclidean 
lattice Lagrangian is the following 

*E( P ,.-.)4e(v.) 



9 2 / A lo / 9 

g 2 aj ^ g 2 a xl<J 



l - J2 Tr(U ij (x) + Uj j (x)-2). (3) 



9 2(X x^<j 
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Here Uij(x, t) denotes a space-like plaquette where the first link goes in direc- 
tion i and the second link goes in direction j, and Uio(x, t) denotes the corre- 
sponding time-like plaquette. One should note that P i0 = P 0i and = Pji 
for plaquettes. In the temporal gauge U (x,t) = 1, the time-like plaquette 
becomes a function only of the link variables Ui(x, t) (i = 1, 2, 3) and we have 

P l0 = ^Re Tr{U l0 ) = ^Re Tr (U t (x,t)Uj (x,t + a t )) . (4) 

We want to construct a classical Lagrangian defining trajectories of gener- 
alized coordinates Ui(x,t) and generalized velocities, where the variable t is 
now continuous. We assume that the action corresponding to this Lagrangian 
is given by the continuum limit a t — > of the lattice action which also yields 
the dependence on the generalized velocities. 

In order to construct this Lagrangian, we introduce a Taylor expansion 
in time and write for a fixed (i, x) up to errors of 0(a^) 

Pio-1 = ^-ReTr[U(t)U\t + a t )-l] 

= ^fe Tr (u{t)\U\t) + a t U\t) + £fr(t)] - l) 
a? „ „ r „, . ''4- / sn a 2 



Re Tr[U(t)W(t)} = --l-Tr[U\t)U(t)} 



2N C 2N C 
= -^TrW)m\- (5) 

Here, we denote U(t) = Ui(x) = Ui(x,t), and we have introduced a gener- 
alized velocity q (corresponding to the angular velocities of the classical top 
theory) which is an element of the SU(N C ) Lie algebra, 



q a (x)\ a = -iUi(x)U}(x) = iUi{x)U}{x). (6) 



The SU(N C ) generators A a are normalized to tr(X a X b ) = 5 ab /2. Going to the 
limit at — > and performing a transition from Euclidean to Minkowski space 
(t — > it) yields the classical lattice Lagrangian 

E + 4-E Tr(Uv + 4 - 2). (7) 



JM 2o 2 ^ v ~ / " 11 v ~ / ' fl 2 fl 
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Here we denote Uij = Uij(x) = Uij(x, t). For a classical canonical formulation 
we introduce the conjugate momenta 

= B«(x)A« = (8) 

The standard Legendre transformation leads then to the following classical 
lattice Hamiltonion 



<9gf (x) 



= f n Y,E?(x)E?(x)-^- 52^(^ + 4-2). (9) 

za x,i & a x,i<j 

Recalling —>■ gaA, we convince ourselves that Ei « a 2 Ai/g. Therefore 
Ei(x) is the approximated color-electric field on the lattice. 

To quantize this classical theory, we proceed according to the rules of 
quantization of the classical top theory ||14| . This results in the prescription 
that the quantum mechanical states are functions of the link variables Ui(x) 
and that the canonical conjugates Ei(x) become the operators of infinitesimal 
left multiplication. A generalization of the standard quantum mechanical 
formula x + a = Q- l P a xe l ^ a yields for the link variables 

e^ xa U 3 {x) = e-^ E ?^U 3 {x)e^ E ^\ (10) 

Variables corresponding to different lattice links are considered to be inde- 
pendent. This yields the commutation relations 

[U t (x),E«(y)} = X^x^ySij, 

[U}(x),E«(y)) = -U}(x)X a 5 x , y 5 l , r (11) 

Since the operators e lE i ^ A yield a representation of the gauge group 
SU(N C ), we obtain for the E 3 (x) the commutation relations of the Lie al- 
gebra 

[E«(x),Ef(x)]=ir^Ej(x). (12) 

The quantization of the classical Hamiltonian, Eq.(^), by use of the commu- 
tation relations Eq. flTTl) gives the standard quantum Hamiltonian of Kogut 
and Susskind. 
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2.2 The transfer matrix method 



The construction of the Kogut Susskind Hamiltonian from the Wilson action 
via the transfer matrix method has been first established by Creutz|l3| (see 
also Ref.[n]). Here we recall the basic steps which may be used also for the 
construction of the improved Hamiltonian as discussed below. 

S E = [dtL E = J2 at[L (q(t), q(t + a t )) + Li(g(t))] + 0(a 2 t ). (13) 

Hereby, L Q is kinetic part of the Lagrangian which couples the system at time 
t to that at time t + a t . Invoking the Baker-Cambell-Hausdorf formula and 
going to the limit a t —>■ ||12|| , the Hamiltonian is eventually given by 

H = H + L U (14) 

where the nontrivial part Ho is related to Lq via the functional integral kernel 
of the corresponding time evolution operator (transfer matrix formalism). 

Let us recall that relation for the simple example of standard one-body 
quantum mechanics of free motion |12| where one has (we put the mass 
m = 1) 

L (q',q) = ^(q' - q)\ (15) 

with q' = q(t + at), q = q{t)- The discrete time-evolution, which relates the 
generator Hq to the kernel Lq, is given by 



c 



-a t H 



^){q) = N [ dq'e- atLo ^ q,) ^(q'), (16) 



where N is some unimportant normalization factor. Using 

Hq')=e {q '- q)v Hq), (17) 

yields 

e -atH = N J dq > e -atL (q,q>) e (<?'-<?)V_ ^ 

In this case, and also for the case of scalar field theory [12], this integral is 
analytically computable for finite a t . It yields the usual result H = —A/2. 
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In order to clarify the notations and the particularities for non-abelian 
gauge field theories, we first recall how to obtain from the transfer matrix 
method the standard Kogut Susskind Hamiltonian. We start by the decom- 
position of the action as given by Eq.([l]) and introduce the temporal gauge. 
The corresponding Lagrangian is given by 



L {U{t + a t ),U{t)) = ^J2 Tr ( 2 - U i( x ' t + a t) U i( x ' t )- h - c -)' 

Li(U(t)) = 4- E Tr{2-U ij (x,t)-h.c). (19) 

9 a x,i<3 

Here L corresponds to the kinetic part and L\ to the potential part, respec- 
tively. The kinetic part of the Lagrangian is given by the plaquettes involving 
different times. Using the notation analogous to Eq. ([T5|) , 



a 

2 . 

V*j = ULjVL, (20) 



9 a t 



where U corresponds to the time slice t and U' to the time slice t + a t , 
respectively. It is well known that the quantum mechanics of SU(N C ) gauge 



theory and that of the quantum mechanical top are closely related 0. |T4 
Thus the relation between the wave function at time slice t and time slice 
t + at, in analogy to Eq. (|TTD , involves the standard color electric operators 



E a ,a = l,..,Ni-l, 

V(U') = e^ aEa m(U). (21) 

where the numbers u a are the parameters of a group element gy £ SU(N C ) 
such that 

gv = e iuJ " xa = U'U- 1 = V. (22) 



The analogy to Eq.QTS]) the Hamiltonian H in the case of one link variable 
V = U'W is given by 

e -«tH =N J dU , e -<HL (V) e i^E^ (23) 



where L (V) is given by Eq. (|20|) . We use the invariance of the Haar measure 
yielding dU' = dV and note that like in Eqs.( |TT| , |2~T|) , the operators E a have to 
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be treated as commuting with U and U' . The integral in Eq.(p3[) can not be 
evaluated analytically for finite time translations at- However, for the deter- 
mination of H one may consider the limit a t — > 0. In this case, the variables 
V approach the identity and it is legitimate to use the approximation for the 
L term 



Tr{V + V ] ) = 2Tr(cosAV> c 



N c - \j 2 + 0(lo 4 ) 



(24) 



Writing the group integral J dU as Haar measure / Y\ a du a detj ac yields for 
Eq . fl23|) a Gaussian integral in analogy to Eq . (|T8|) , with the exponent 

a -uo a uo a + iE a u a = -^(u; a + iE a< ^) 2 + ^E a E a . (25) 



2a t g Zcbtg a 2a 

Taking the sum over the space-like links this reproduces the standard kinetic 
term of the Kogut Susskind Hamiltonian 

^o = ^EWW- ( 26 ) 

3 Classical improvement of Kogut Susskind 
Hamiltonian 

3.1 Continuum behavior of classical improved action 

The Wilson action reproduces the classical continuum action only up to errors 
of 0(a 2 ). It is possible to add to the Wilson action new terms such that these 
0(a 2 ) errors are canceled [VJ\, [TB], M. In order to construct the corresponding 



improved Hamiltonian, one needs a generalization to lattices with at ^ a s = 
a. We first discuss the classical continuum behavior of the Wilson action. 
For space-like plaquettes one has (see Ref.|l|) 



Pa = Y Re Tr{u ^ 1 ~ Cv" [ Tr Fv F d + y^\F ld (v 2 + v 2 )F l3 ] 



(27) 



For convenience, the continuum limit of a loop is expressed in terms of the 
field strength tensor and its covariant derivative the center of the loop Xq. 
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According to Ref.jl]] for time-like plaquettes one has to consider the path 
ordered integral 



r ra/2 rat/2 1 

fA-dx -> / dxi / (it[F i0 (a;o) + -xvXjiiVyV^F^ 

J J-a/2 J-at/2 2 



ra/2 j-at/2 

~'iO)\x=xo\ 

'-a/2 J— a t /" 

3 3 

aa t F i0 {x ) + ^(V^)F t0 (x ) + ^(V 2 )F l0 (x ). (28) 



The last term can be neglected since a t « a. Therefore 

P t0 = !-ReTr{U l0 )^^ReTr[l- l -{gjA-dxf 



v c 



1 Re Tr[l - 9 -A^(F t0 + ±-a 2 V*F i0 f] 



N c 2 y 24 

-,2„2„2 / „2 



1 " ^ f Tr[F l0 F l0 ] + ^Tr{F l0 V*F l0 }\ + 0(a 2 ^). (29) 



In order to compensate these 0(a 2 ) errors, Lepage |I] has proposed to add 
new terms to the Wilson action. One of these terms is given by a rectangular 
loop, 



Rp, = j^-ReTr 



v 
\ 



X 



(30) 



For a space-like loop one has in particular 
1 



Ri 



N, 



Re Tr[Ui(x)Ui(x + ai)Uj(x + 2ai) 



x U}(x + aj + ai)Uj(x + aj)uUx)] 



1 



2 N r . 



4Tr[F tJ F l3 ] + -Tr[F^V\ + V])F l3 ] 



(31) 



Considering time-like loops, there are two possibilities. First, one has a 2a x at 
rectangular loop 



R i0 = —ReTr 




t 



(32) 



x 
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yielding 



RiO 



X 



— — Re Tr[Ui(x, t)Ui(x + ai, t)Uo(x + 2ai, t) 
Uj (x + ai,t + a t )Uj(x, t + a t )Ul{x, t)\ 




Secondly, one has a 2a t x a rectangular loop 





(34) 



This term corresponds to advancing two steps in time direction. The 
conventional transfer matrix corresponds to an advance of a single step in 
time direction. Thus it is not compatible with the definition of the transfer 
matrix. We may disregard this term because the improvement terms in the 
Lagrangian are not uniquely determined 0]. Taking into account only the 
first term is sufficient. 

Therefore, we make the following ansatz for the classically improved Eu- 
clidean lattice Lagrangian Q. 



2N c a ^ , P i0 + P 0i , , 

L t = YTzA c i o h C 2J Ri J + const. 

9 a t x>i 1 

C[ = 4/3 and C 2 = -1/12, 
L s = J2[C l ^±^ + C 2 (R tJ +R Jl )}+ const. 

Ci = 5/3 and C 2 = -1/12, 
L E = L t + L s . (35) 
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3.2 Improved Hamiltonian via Legendre transforma- 
tion 



Now we proceed as in sect. fl2.lp to construct a classical Lagrange function 
in Minkowski space in terms of the generalized coordinates Ui(x) and the 
generalized velocities qi(x) as defined in Eq.(|B|). Working in the temporal 
gauge and denoting U(t) = Ui(x + ai,t) yields 

R i0 = ^ReTr[U(t)W(t + a t )U"<(t + a t )U(t)] 
-> j^-Re Tr{U{t)[U\t) + a t lj\t) + ^u\t)\ 
x [u\t) + a t U\t) + ^W(t)]U(t)} 

i 2 2 ± 

= W Re Tr[°±U{t)fj\t) + ^U\t)U{t) + alU (t)(j\t)U\t)U (t)} 
+ const. 

—> -^-a^Tr[]-qi(x + ai)qi(x + ai) + ^qi(x)qi(x) + Qi{x)qi(x + ai)] 

+ const., (36) 

where we have introduced the variable 

Qi{x) = U t (xyq t (x)U t (x). (37) 

This gives the following classical improved lattice Lagrangian in Minkowski 
space 

L M = 4 E Tr K C i + 2C' 2 ) qi (x) qi (x) + 2C' 2 Q i (x)q i (x + ai)} - L s . (38) 

This Lagrangian can be written in the form 

1 a 



m = ~^J2q a (t)M ap (U(t))q p (t) - L 8 , 
where a=(x,i,a), p= (y,j,P), 



+ 4C' 2 5(x,y -ai)5 iJ Tr[U}(x)\ a U i (x)\> 3 }. (39) 



rt, 
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The matrix M is not symmetric. However, it can be shown that only the 
symmetric part of M will contribute to the Hamiltonian. Thus we introduce 

L%> = 6(x,y-al)5 t ,Tr[U^ y \ a U x ^} 

+ 5(x,y + al)5 h3 Tr[U^ x \PU y ^ x \% (40) 



which allows to write 

1 
2 



MZ m = tt( m + m %p = (C[ + 2C^ CT , p + 2C^L^. (41) 



Inspection shows that L S3/m and hence M sym are real, symmetric matrices. 
Then the Lagrangian reads 

L m = ~E?,W M^{U(t)) q p (t) - L s . (42) 
Via Legendre transformation, the classical improved Hamiltonian is obtained 

#° = ~E^(M"»™)-X 

2iV c r Pi - + Pi 
V = — ^ Ci 13 — — + C 2 {Rij + Rji 
9 2 a X %1 2 



(43) 



The color-electric field E a is given by the conjugate momentum, being related 
to the generalized velocity q a via, 

^g^M"V (44, 

The color-electric field, obeys commutation relations with the link variables 
given by Eq. ([TTj) . 



3.2.1 Hopping expansion and algebraic properties of M t 



sym 



Taking a closer look to the kinetic part of the improved Hamiltonian reveals 
that via M~ ym an infinite number of terms enters into the Hamiltonian. In 
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analogy to the hopping parameter expansion JTT] , which expresses the prop- 
agator in terms of powers of a hopping matrix, we introduce if sym 



M. 



sym 

m 



m [1 + K sym ] = m [1 + k L 
7 

C 1 + 2C 2 = — , 



sym] 



k, 



2C' 



o 



1 

'V 



to obtain 



M" 1 

sym 



1 

m 



1 — if + if — if d + ■ ■ 

x ly sym i ly sym ly sym ~ 



(45) 
(46) 



While if sym involves only link variables between next neighbor lattice sites, 
higher powers of K sym involve links extending over several lattice sites. Using 
the notation 



U x ^ x+(N+1)a - = Ui(x)Ui(x + ai)---Ui(x + Nai), 
we generalize the definition of L sym to 

L (N)sym = §( x ^y _ Nal)5ijTr[Ul_ ty \ a U x ^y\ (S ] 
+ 8(x, y + Nat)5 hJ Tr[U^ x \PU y ^ x \% 
N =0,1,2,-.., 



(47) 



(48) 



where LfJ m = 1, and Lg m = L 

obeys the following product rule 



„, y „, .... ...... -sym- A little algebra shows that the matrix 



j (p) r (?) 

sym sym 



Ir>+?) + ir/b-gl) 

2 sym ' 2 sym 



(49) 



Thus we obtain for the lowest powers of if, 



Ksym 


— k L^J m 


5 


if 2 


- h 2 

— K, Q 


2 


(0) , ir (2) 

sj/m ~ 2 sym 


if 3 


— ^0 


[<§ 


+ ~)LM + 

1 ^ / sj/m 1 


if 4 

sym 


— tVQ 




"I - q )^- / sj/r?i "I - 



^sym 
2 sym 



sym 
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K 



K 



sym 



+ 



K 



k 6 
K 



f i + L\Tp) + 4 + — u,< 3 ) + — rS^ 

^2 g ' sym ^ 4 10 sym sym 

4 + 4)l(°) + (I + 1 + 4 + 4 )L (2) 



sj/ra 



+ 4 + 4 )L (4) + 

8 \Q J s v m 



32 S2/m 



It has the general structure 



(50) 



t,™ \^ K-( n ) r (p) 

^0 p ^sym- 



p=0 



(51) 



The coefficients of lowest order are 



jo) 

K (l) 

(2) 



(1) 



o, 4 

I K (2) 

2' 1 



o, 4 2 > 



1 

2' 



(52) 



The coefficients tzj^ vanish except when n and p are both even or both odd. 
Using Eq. fl5T|) , we express M~y m by 



M" 1 

sym 



■y oo 

^P-L&yr 



m 



symi 



p=0 



»v = E(-^o) n 4 n) - 



(53) 



n=p 



As result, starting from an improved Lagrangian with a finite number of 
terms, one obtains for the improved Hamiltonian an expression given by an 
infinite number of terms. 

In the following we will explore more of the algebraic structure of M sym 
and obtain analytic expressions for the hopping expansion coefficients 
This will be useful in what follows. We introduce 



J CT , P = 25{x,y- ai)8 i;j Tr[tf \ a U x ^ y \P\. 



(54) 
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A little algebra shows that 



Jl P = 25(x,y- nad)5 id Tr[UU y X a U x ^], 



(55) 



and 



J J* = J 1 J = l, 



(56) 



i.e., J is a real, orthogonal matrix. Comparison with Eqs. (|4C| , j48|) shows 



J- 1 sym 

L (P) 
sym 



K 



sym 



M, 



sym 



m 



1 + |(J+J- 



Using Eq. (|56"D , M sym can be factorized, 



m 



L sym ]_ _|_ £2 

if C is chosen as solution of 

k = 



;i + cj)(i + cj*), 



2C 



l + C 2 ' 



(57) 



(58) 



(59) 



Solutions are C = -7 ± 4V3. Note that J being a real, orthogonal matrix, 
which has eigenvalues of modulus one, and \C\ ^ 1, thus the matrix M sym 
can be inverted and M~y m is well defined. Moreover, we note that M sym is 



a positive definite matrix. This can be seen directly from Eq. (|58|) , which 
factorizes M sym into a matrix times its Hermitian conjugate. Also, a lower 
bound can be estimated using Eq. (|57|) . J being orthogonal implies || J\\ = 1. 
Thus Rj defined by Rj — |(J + </'), being a real, symmetric matrix like 
obeys | \Rj\ \ < 1. Then an arbitrary state vector of unit norm yields 



symi 



< 



\R:, 



> | < 1. Then Eq . (|5^) implies 



< 6\M. 



sym | 



> = m + m k < (j)\Rj\(j) > 

= I - I < <P\Rj\<P > >i, 



(60) 
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showing also that M sym is positive. To summarize the properties of M sym , 
this is a real, symmetric, positive definite and non-singular matrix. This 
property is needed for the construction of the Hamiltonian via the transfer- 
matrix, in particular for doing the Gaussian integral. 

Factorization of M sym , via Eq. fl58|) , allows to express the kinetic energy 
term H of the Hamiltonian, Eq.(^), as follows 



Note that this is an expansion in terms of C and J. 

Analytic expressions for the coefficients of the hopping expansion can be 
obtained in the following way: 



H 



= - Tr Y [EAx) - CUi{x)Ei{x + ai)Uj(x) 

a m o tt 

+ C 2 Ui(x)Ui(x + ai)Ei(x + 2ai)uj(x + ai)uj(x) - 




2 



(61) 



K: 



sym 



71 




(62) 



Because this expression is a symmetric matrix and making use of Eqs. fl56| , |57| ), 



one may write 



t, n u n 

Km = (i) n E ^(J p + (J'T) = O n E «i n) 2^i- (63) 




Comparison of coefficients yields 



p = : a 



P 




(64) 
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Comparison with Eq.([H]) eventually yields for the hopping expansion coeffi- 
cients k^ 1 ' the following expression, 



K {n) = 2 -n+l a (»). 



(65) 



3.3 Improved Hamiltonian via transfer matrix 

We start from the classically improved Euclidean Lagrangian, given by 
Eq. (|35|) . It is built from space-like plaquettes Py, time-like plaquettes Pio 
and corresponding rectangular loops Rij and Riq. We now want to show 
that the transfer matrix method yields the same Hamiltonian as has been 
obtained in the previous section via Legendre transformation. Let us con- 
sider the time-like part of the Lagrangian, which yields the kinetic part of 
the Hamiltonian. The space-like part yields the potential part in a trivial 
way. Using the temporal gauge, one has 

P 1 



2N t 
1 

2N r 



■Tr 
■Tr 



Ui(x, t)U}(x, t + at) + Ui(x, t + a t )Jj}{x, t) 

Vi{x,t) + v}(x,tj 



(66) 



using the notation Vi(x, t) = Ui(x,t + a t )Uj(x,t). Similarly one obtains for 
the rectangular loop 

1 



R 



iO 



2N r 



-Tr 



Vi(x, t)Ui(x, t)Vi(x + ai, t)Uj(x, t) 



+ Ui{x, t)V?(x + ai, t)U}(x, t)V?(x, t) 



(67) 



The Hamiltonian is defined via the transfer matrix like in Eq. (p3|) . Because 
we consider a t — > 0, the group integral will be dominated by group elements 
of SU(N C ) in the neighborhood of the unit element. Thus one can expand 
the group elements Vi(x,t) in a Taylor series of the Lie group parameters 



Vi(x,t) = exp[iu X)i (t)] = 1 + iu X)i (t) - \j x ,i(t) 2 + 0(u 3 ), 
where we denote u> Xti (t) = J2 a ^ a ■ Thus we arrive at 



(68) 



Tr\V t (x,t) + V?(x,t)] = 2N C - - Y.^M + °^ 



(69) 
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and hence in the notation of Eg. (|39|) , 

Y,Tr[Vi{x,t) + V?(x,t)] = -~5> ff (*) 5 °p u pV) + °(^ 3 ) + const - ( 7 °) 



ap 



Carrying out the corresponding steps for the rectangular term, one obtains 
Y^Tr \Vi(x, t)Ui(x, t)Vi(x + at), t)U\{x, t) 



+ Ui(x, t)Vi(x + ai), t)Uj(x,t)V^(x,t) 



ap 



$op 2 ^ ap 



uJi) + 0(u 3 ) + const, 



(71) 



where J is given by Eq . (|54| ) . Inserting this into the time-like Lagrangian Lt 
yields eventually 

1 " - \ ,..>_tt) K ". .... •>( •'.)()_ ... if ",., T(j 



Y ^{t) \{c[ + 2C' 2 )5 ap + AC 2 ry p m } Up (t). 



(72) 



One should note that the matrix J is not symmetric. However, to the La- 
grangian only the symmetric part J sym — (J+J t )/2 contributes. Note further 
that J sym = L sym and M sym = (C[ + 2C' 2 )1 +4C' 2 J sym , being real, symmetric 
matrices. Thus we arrive at 

£. = ^EMt)W)- ( 73 ) 



2 ah 2 



ap 



The transfer matrix is then given by 
exp[-ai/7 + 0{a 2 t )\ 
lUdViix)} eM-aM^x))} exp[* Y (*)] 



= iVexp 
Thus we obtain 



Y[ du a det Jac \ exp 
1 



1 a 



2a t g 2 



Y^M^u. + iY^E, 



ap 



YM— 2 M s n^Ep 

a t g 



ap 



H = ~^lY E v( M sym)apE p . 



(74) 



(75) 



ap 



in agreement with the result, Eq. 



obtained via Legendre transformation. 
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4 Improved Hamiltonian given by finite num- 
ber of terms 

As was shown in the previous section, the kinetic energy of the classical im- 
proved Hamiltonian obtained directly from Lepage's action is given by an 
infinitive series of terms. Even though the series is rapidly convergent, such 
an Hamiltonian is too complicated for a practical calculation. Recalling that 
the purpose of classical improvement is to push the 0(a 2 ) error to order 
0(a 4 ), we show here how to construct a simpler improved Hamiltonian cor- 
responding to a finite number of terms to achieve such a goal. In the previous 
section we have seen that the infinite number of terms in the Hamiltonian 
arises due to the inversion of the matrix M sym , which itself has only a finite 
number of terms. Thus it is plausible that in order to obtain a Hamiltonian 
given by a finite number of terms, one needs to start from a Lagrangian 
corresponding to a matrix M sym with an infinite number of terms. Such a 
construction is possible, because the Lagrangian leading to improvement is 
not unique. We start by considering the following type of Wilson loop, which 
emerges as a generalization of the 2a x a t loop to a (n + l)a x a t loop 
parallel transporter R n i t o given by 

| (76) 

x + nai 

In the temporal gauge it corresponds to the expression 

Rni,o = —Re Tr[Ui(x)Ui(x + ai) ■ ■ ■ Ui(x + (n — l)ai)Ui(x + nai) 

x U}(x + nai, t + a t )U}(x + (n - l)ai) ■ ■ ■ U}(x + ai)U-(x,t + a t )]. 

(77) 

Note that for n = 1, R ni:0 coincides with R i0 . The path-ordered integral of 
such a Wilson loop is given by 

j>A-dx -> J dxidt{F i0 (x ) + -XitXj'iVitVj, 

Fio) \ x=xq\ 



Rnifi — 



-Re Tr 



x 
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1 



(n+l)a/2 



(n-l)o/2 



2aa t F i0 (a;o) + ^a t V^F i0 (x )[ 
3 n 2 _)_ i 

2aa t F i0 (x ) H — — a t a 3 T> 2 F i0 (x ). 



oc di^jc I 



-(n-l)a/2 
(n+l)a/2 



12 



x 2 fix] 



(78) 



Therefore, we obtain the following continuum behavior for the above Wilson 
loop 



Fn.i 



nifi 



±. Re Tr[l - ^(f A ■ dxf 



2„2„2 



g a a t 
2N r 



(4Tr[F iQ F i0 \ + (n 2 + l/3)a 2 Tr[F i0 V 2 F 



(79) 



One verifies for n = 1 that Eq . (|79|) coincides with Eq . fl33|) , as should be. 
We make the following ansatz for the Euclidean lattice Lagrangian 



2N c a A , 



9 2 a 2 



n=l 



(80) 



In order that the usual continuum limit of the Lagrangian is obtained and 
the 0(a 2 ) error is canceled, we imply from the continuum behavior of Pm, 
Eq . (|29|) , and of R n i,o, Eq . (|79"D , that the following conditions hold, 



A'[B' + 4 J2 C' n ) = 1, 

n=l 

r>l oo i 



(81) 

" !,i -0. (82) 
We have deliberately introduced the coefficient B'. Choosing 

oo 

B' = 1 - 2 ]T C' n , (83) 



n=l 



results in a simple expression of the Lagrangian expressed in terms of gener- 
alized coordinates and velocities. Using 



nC - c dv [C do ) t 1 ~ (i - cy ' 



(84) 



n=l 
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we obtain 

i -a 



A' = 



1 + C' 



B> = (85) 

and C" is a root of 

C ,3 + 11C" 2 + 11C+ 1 = 0. (86) 
This equation has three real roots, given by 

C = -1, 

C" = -5±2>/6. (87) 

The root closest to zero is C' = -5+2^ = -0.101021 • • In order to obtain 
the kinetic energy, we express R n i,o in terms of generalized coordinates and 
velocities, 



—Re Tr 



[Uj(x) + a t U}{x) + ^a 2 t Uj(x)]U t (x)U t (x + ai) 



x Ui(x + (n — l)ai)Ui(x + nai)[ll}(x + nai) + a t Uj(x + nai) 
+ -a 2 t)-(x + nai)]l/J(x + (n — l)ai) ■ ■ ■ U}{x + ai) 
11 1 - 

— a\Qi(x + (n — l)ai)(ji(x + nai)], 

(88) 



where we have introduced 
Qi(x + (n — l)ai) 



= Uj(x + (n- l)ai) ■■■Uj(x + ai)U}(x)qi(x)Ui(x) 
x Ui(x + ai) ■ ■ ■ Ui(x + (n — l)ai). 

(89) 
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Note again, for n = 1, Qi(x + (n — l)az) coincides with Qi(- C ) defined in 
Eq. (]37f) . Thus we can write the time- like part of the Minkowski lattice 
Lagrangian, 



- M 



qi{x)qi{x) + 2 ]T C' n Q % {x + (n - 1) cm) ft (a: + 



new 



n=l 



.(90) 



To find the kinetic energy of the improved Hamiltonian, it is convenient to 
express this in terms of the matrix J, defined by Eq . fl54f) , 



where 



M' = A' 



2g 2 cj p 



1 + 2 ]T C /n J n 

n=l 



.4 



,1 + C"J 

1-C'J' 



The color-electric field is expressed as 

<9L a 



E a 



, 2 E M TV 



(91) 



(92) 



(93) 



where 



ri-c 



'A 2 



i i _ r" 2 

2 y ' (1-C"J)(1-C"J*) (1 + C" 2 ) - C"(J + J*) 

(94) 

If we choose C such that |C"| ^ 1, e.g. C£ = -0.101021 • ■ ■, then M' sym is 
a real, symmetric, positive and non-singular matrix. Finally, we obtain the 
corresponding kinetic energy of the improved Hamiltonian, given by 



2a 



2a 

g" 



a.p 



E 

c,p L 

E 

a,p 



— Tr V 
a ti 



1 + C' 2 
[1 - C') 2 

1 + C' 2 

{i-cy 

1 + C' 2 



[l-C 
2C 



2C 



—E a (J + J l ) atP E p 



(l-C 



E a 5 at pEp — 
E a 5 ai pEp — 

Ei(x)Ei(x) - —p=^ Ui(x) ] 'Ei(x)Ui(x)Ei(x + ai) 
(1 — O J 



7\2 



(95) 
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It consists of only two terms, which makes it convenient for practical calcu- 
lations. 



5 Tadpole Improvement 

In the preceding section, we derived a classically improved Hamiltonian for 
gluons with 0(a 2 ) corrections. A very important step of the improvement 
program is to take into account quantum corrections corresponding to tad- 
pole diagrams. Without such improvement, only part of the 0(a 2 ) errors are 
canceled. According to Lepage and Mackenzie, most of the tadpole contri- 
butions can be removed by dividing each link operator by the mean 
of the link. For asymmetric lattices, a t << a s , and small enough a t we have 
u t = 1 for time-like directions. In the Hamiltonian formulation, the mean u s 
of a space-like link is defined by 



Us 



<n|iW 1/4 > (96) 



where \Q) is the vacuum of the improved Hamiltonian. Thus tadpole im- 
provement of the lattice Lagrangian L = L t + L s , where L t is given by 
Eq.flgQP and L s by Eq.(R3), corresponds to the replacements 



Pij — * Pijl U i 
Rij Rij I XL 



it 

Rnifi Rnifl/U 2 , n+ . (97) 



This is equivalent to the following replacement of constants 



d -> d/u 
C 2 -> C 2 /u 



9t -> 9tu s , 

.2 



C - C'/ul (98) 

where we put g = g t in Eq.(^). For the transition to the Hamiltonian these 
redefinitions of the coefficients can be taken over yielding for the "two-term" 
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tadpole improved Hamiltonian (C = 



-0.101021) 



H 

H 



V 



H + V, 



2„,2 



9t U 



x,i 

2C'/u 2 s 
(1 - C'/ul) 

s*" X,i<j . 



1 + C ' 2/u ' : E,(x)E i (x) 



9 z s a 



[i-c/uir 



- Ui(x) ] ' Ei(x)Ui(x)Ei(x + ai) 

C x Pj j + Pjj c 2( , 



U 



u 



(99) 



Here, we have introduced different couplings in the kinetic and potential 
terms in order to allow for a "speed of light" correction as discussed in Ref . [^3| 
(see below). 



6 Further perturbative improvement 

Tadpoles have been identified as an essential part of the problem when ap- 
proaching the continuum limit of quantum field theory on the lattice. A 
systematic perturbative calculation on the lattice has been performed by 
Liischer and Weiszp0[. This leads to the determination of additional terms in 



the Lagrangian needed to compensate errors. It turns out that such a further 
improved Lagrangian (for details see Ref. |0| ) contains the same plaquettes 
and planar rectangle loop terms which occurred before, but with suitably 
redefined coefficients, plus a new term, being a non-planar "parallelogram" 
loop, given by 



C nun — -—-Ra Tr 



c 



(100) 



It corresponds to 



C x ^ ua — ^ Re Tr(U x ^U x + a ^ tV U x ^ a ^j tau ^U x+a +a ^^U x+a ^ v U xa ^ (101) 



The corresponding term occurring in the Lagrangian is proportional to 

E (102) 

x,)i<v<a 
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The structure of the corresponding improved Hamiltonian can be inferred 
from the improved Lagrangian as before: One introduces different lattice 
spacings a s = a and a t and constructs the Hamiltonian by Legendre transfor- 
mation and canonical quantization. Here, we refrain from discussing details 
and only give the general structure of emerging Hamiltonian. 

(1) The plaquette and planar rectangle loop terms will give a part of the im- 
proved Hamiltonian which has the same form as before, only the weights of 
the individual terms will be different. 

(2) The space-like parallelogram loop terms (i.e. \ivo~ space-like) will yield a 
corresponding additional term in the potential part of the Hamiltonian. 

(3) The time-like parallelogram loop terms (where either \x or v or a is time- 
like, the other two indices being space-like) produce a large number of differ- 
ent contributions to the Hamiltonian (with well defined weights). The final 
result for the improved Lagrangian has the structure 

L = L t (q,U)+L s (U), (103) 

with 



T - a 



(104) 



q*M ap (U)q p + ]T A a (U)q a + h.c. 

. &P a 

A new feature is the occurrence of a term linear in q. As before, M(U) is a 
symmetric matrix of the form 

M = l + M, (105) 

allowing the definition of M _1 by a geometric series expansion. Legendre 
transformation and quantization yields a Hamiltonian of the structure 

H = H + V, (106) 

with 

Ho = fj2 {EoKp'Ep - {A a M^E p + h.c.) - 2A a M^A p ] . (107) 



ap 



7 Discussion 

For the purpose of a numerical calculation, in particular for a comparison 
with lattice Monte Carlo results, the following points are important: 
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(1) As discussed in Refs. |22|, f2"3f , the scales related to the regularization of 



the gauge field theory in the Hamiltonian formulation as opposed to the 
Euclidean path integral formulation are different. This difference can be 
accounted for by introducing spacelike (g s ) and time-like couplings (g t ) which 
have a well defined relation to the "Lagrangian coupling" g. In one-loop 
approximation this relation is of the type 

4 = ^+ c /- ( ios ) 

where c„ depends on the space-time dimension and on the type of the gauge 



group and is given in detail in Refs. [p2| , [24| , £3 



(2) Because of this difference in the nature of the lattice regularization, all 
perturbative calculations which determine some non-classical improvement 
in the sense of Luscher-Weisz have to be redone. Such a calculation can be 
done on an asymmetric Euclidean lattice with a t « a s (see Ref . [Ef ) . 



(3) Tadpole improvement which has been considered by Lepage flTUf in the 



Lagrangian framework corresponds in the Hamiltonian framework to an ex- 
pression given by Eq. fl99|) . 

(4) A systematic determination of the Luscher-Weisz improvement terms on 
asymmetric lattices in the Hamiltonian framework has still to be done. Since 
these additional corrections turn out to be small in the standard Euclidean 
framework (see Ref.Jl9|) - the most important correction coming from the 



inclusion of the tadpole terms - in should be worthwhile to work with the 
improved Hamiltonian given by Eq. (|99|) , e.g., for the numerical simulation of 
glueballs. 

To summarize, we have investigated in this paper two schemes of improve- 
ment of the Kogut Susskind Hamiltonian: If one starts from Lepage's La- 
grangian, which is preferable for Monte Carlo simulations in the Lagrangian 
formulation, the corresponding Hamiltonian is given by an infinite series of 
terms which contain terms with arbitrary long range. In contrast, we have 
shown that by starting from a suitable Lagrangian with an infinite number 
of terms , one can get an improved Hamiltonian consisting of a finite small 
number of terms. This should be preferable for numerical computations in 
the Hamiltonian framework. 
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